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Abstract. By means of the Bessel operator a polynomial sequence is constructed to which 
several properties are given. Among them, its explicit expression, the connection with the 
Euler numbers, its integral representation via the Kontorovich-Lebedev transform. Despite its 
non-orthogonality, it is possible to associate to the canonical element of its dual sequence a 
positive-definite measure as long as certain stronger constraints are imposed. 



1. Introduction and preliminaries 

The modified Bessel function of third kind (also known as the Macdonald function, specially 
in russian literature) Kj X (x) of the argument x > and the pure imaginary subscript i% is an 
eigenfunction of the (Bessel) operator 

. 2 a ] d 2 d 2 d 2 /- -,\ 

=xr'—x—x— = —x —^—x——+x (1-1) 
dx dx dx 1 dx 

for the associated eigenvalues T 2 , i.e., 

sJK H (x) = T 1 K H {x) (1.2) 
and naturally providing the identity on the integral powers of the Bessel operator srf 

srf n K H {x) = T 2 "K iT (x) , n € N, (1.3) 

inductively defined. Throughout this paper, N will denote the set of all positive integers, 
No = NU{0}, whereas R and C will denote, respectively, the field of the real and complex 
numbers. By R + and R ( | we respectively mean the set of all positive and nonnegative real 
numbers. Further notations are introduced as needed throughout the text. 

Defined also by the cosine Fourier transform 

/■ CO 

K H {x) = / e - AXOsh(l,) cos(ru)du , x e R+, iel 
Jo 
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the function K\ x is real valued and represents a kernel of the following operator of the Kontorovich- 
Lebedev transformation |19| . given by the formula 

fOO 

Jtlxif} = / K iz {x)f{x)dx , (1.4) 
Jo 

is an isometric isomorphism (see [17]) between the two Hilbert spaces L2(M. + ;xdx) and L2(M + ;Tsinh 
Moreover the inverse formula holds 

f( x ) = — zsmhnz K ix {x)J^ % [f\ dy (1.5) 

7l z X Jo 

where the convergence of the latter integral is understood with respect to the norm in L2(^ + ;xdx). 

The Kontorovich-Lebedev transform has been used in many applications including, for instance, 
fluid mechanics, quantum and nano-optics and plasmonics. This transform, where integration 
occurs over the index of the function rather than over the argument, proves to be useful in 
solving the resulting differential equations when modeling optical or electronic response of such 
problems. 

Recently [16J, the action of any power of the operator over e~ x was investigated. Precisely, 
the sequence {e.~ x .saf n e x } n ^ was shown to be a polynomial sequence (PS), spanning the vector 
space - the set of all polynomials with coefficients in C - since each of its elements has 
exactly degree n. However, this sequence is orthogonal with respect to the singular measure 5 
(the Dirac measure). In the present work, we propose instead to investigate the functions 

e- x x- a £/ n e x x a , n G N , 

which end up to be a polynomials of exactly degree n as long as 9^(a) > —1/2, as it will be 
revealed in Section [2 In this case the sequence {p n (-',cc)} n ^o defined by 

Pn {x;a) = {-l) n e c x- a £^ n t- x x a , ?2GN , (1.6) 

is indeed a PS. The properties of {p n {-\ o:)} n j»o will be thoroughly revealed Section [2] where 
a generating function, the connection with the generalized Euler polynomials or the Bernoulli 
polynomials (when the parameter a ranges in N) will be expounded. On the grounds of these 
developments lies the Kontorovich-Lebedev transform. The differential-difference equations 
over the polynomials p n {x;a) are crucial for the subsequent developments related to the cor- 
responding dual sequence of monic polynomial sequence (MPS) {P n (x;a) := a~ l p n (x\a)} n ^o 
where a n / are the leading coefficients of the polynomials p n (x;a), n G No- 

The research about the MPS {P n (x; ot)}„^o will then proceed toward the analysis of the corre- 
sponding dual sequence, {u n (a)} n ^Q, whose elements, called as forms (or linear functionals), 
belong to the dual space 3 s ' of and are defined according to 

(u n (a),P k (-;a)) := 8 n>k , n,k > 0, 

where 8 n ^ represents the Kronecker delta function. Here, by (u,f) we mean the action of 
u G over / G 2? , but a special notation is given to the action over the elements of the 
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canonical sequence {x n } n ^o - the moments of u G (u) n := (u,x n ),n ^ 0. Any element u of 
can be written in a series of any dual sequence {\ n } n ^o of a MPS {5„}„^o 

u=J^(u,B n )\ n . (1.7) 

To infer differential equations for the elements of the dual sequence it is important to recall 
that a linear operator T : — > has a transpose 'T : 9*' — > defined by 

( t T(u),f) = (u,T(f)), ue&>',fe&. (1.8) 

For example, for any form u and any polynomial g, let Du = u' and gu be the forms defined as 
usual by (it',/) := —(u,f) , (gu,f) := (u,gf), where D is the differential operator [5|. Thus, 
D on forms is minus the transpose of the differential operator D on polynomials. 

In Section [3] we will come to the conclusion that it is possible to associate to wo(oc) a definite- 
positive measure as long as 91(a) > 0. This implies the regularity of any form proportional to 
wo(a), a concept that is recalled in detail on page[T2J Concomitantly, the existence of a unique 
MPS regularly orthogonal (hereafter MOPS) with respect to u<${a) is ensured, that is, there 
exists a unique {Q n (x;a)} n ^o such that (uo(cc),Q n (x;a)Q m (x;a)} = k n (a)d nM with k n / for 
any n,m€ Nq. 



2. Algebraic and differential properties 

The functions defined by (11.61) are actually polynomials of exactly degree n, as previously 
claimed and below proved, and therefore {p n (x;a)} n ^o forms a PS. Once this is guaranteed we 
will proceed to obtain an explicit expression for these polynomials p n (x;a) in §2. 11 an integral 
representation for such polynomials by means of the Kontorovich-Lebedev transform in §2.21 
and finally, by inverting this integral transform, to derive a relation between the generalized 
Euler polynomials and p„(x;a) in §2.31 

Lemma 2.1. For any complex number a such thatSi(a) > —1/2 the functions p„{x; a) defined 
in (II. 6|) are actually polynomials of exactly degree n fulfilling 

p n+ i(x;a) =x 2 p"(x;a)-x(2x- l-2a)p' n (x;a) - {{2a + l)x- a 2 ^p n (x;a) , n e No, (2.1) 

and also 

(2a + 1) x p n (x;a + l) = -p n+ \{x\a) + a 2 p n {x\a) , n e No. (2.2) 

Moreover, p„(0; a) = a 2 " and p n (x;a) = (— 2)"(a + l/2) n x' 1 + a n -\{x) withdega n -\ ^n — l for 
any n£No (under the convention a_i(x) = 0) and (y)„ representing the Pochhammer symbol: 

n-l 

{y)n = n (y + t) for n G N and (y) = 1. 

T=0 

Proof. According to the definition (II. 6p , we may write 

p n+k (x;a) = (-l) k e x x- a £f k (e- x x a (e x x- a £f n e- x x a )) , n,keN, 
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providing 

p n+ k{x;a) = {-\) k & x x- a ^ k {&- x x a Pn {x;a)) , n,k e N. (2.3) 
On behalf of the property, 

e x x- a £/(e- x x a f(x)) = -* 2 J^/(*) +x{2x- 1 -2a)^f(x) + ((2a + \)x- a 2 )/(x), (2.4) 

that holds for any analytic function /, the particular choice of k = 1 in (|2.3I) furnishes f|2. 1|) . 
whereas the choice of n = 1 with & varying within N leads to the equality (12. 2D . 

Insofar as computing (|1.6I) for n = and n = 1 we respectively extract that po(x;a) = 1 and 
p\{x\a) = — (2a + \)x + a 2 , then on the basis of (12. ip . by finite induction we derive that 
p n (x;a) is a polynomial that has exactly degree n as long as 9^(a) > —1/2. 

Finally, the particular choice of x = in (12. 2 p provides p n+ i(0;a) = a 2n+2 po(0;a), while from 
(j2.ip we deduce that the leading coefficient c nn (a) of p„(x;a) satisfies c ll+ \ J1+ i(a) = —(2n + 
2a + l)c, VJ (a) for any n G N, whence the result. □ 

Directly from (12. 3h . we obtain 

£f k (e- x x a p n (x;a)) = ^ n {^ x x a p k {x;a)) , n,ken. 

From the insertion of (j2.2p into (12. ip we derive 

(2a + l)p„(x;a + 1) = -x p"(x;a) + (2x- l-2a)p' n (x;a) + (l+2a)p„(x;a), n^O. (2.5) 

The polynomial sequence treated in [16] corresponds to the special case {p,!(x;0)}„^o- The 
introduction of this slight modification on the sequence {/?,j(x;0)}„j»o - more precisely, the 
inclusion of x~ a on the left and of x a on the right hand side - has the merit and pertinency of 
guaranteeing that we will be able to deal with regular forms as long as 9?(a) > 0, as it will be 
analyzed on Section [3] Notwithstanding this advantage, the analysis of {p n (x',cc)} n ^o became 
significantly more hard to deal with. 

Before proceeding into this, we list a few elements of {p n (x;a)} n ^o 
P\(x;a)=x(-2a- l) + cr 

p 2 (x;a)=x 2 (4a(a + 2) + 3)+x(-2a(a(2a + 3) + 2)- l) + a 4 
/7 3 (x;a)=-x 3 (2a + l)(2a + 3)(2a + 5)+x 2 (2a + l)(2a + 3)(3a(a + 2) + 5) 
-x(2a + 1) (a 2 + a + 1) (3a(a + 1) + 1) + a 6 

The polynomials p n {x; a) in the variable x are also polynomials in the variable a, but of degree 
2n for each n € N. This can be deduced from (j2.2p . 

2.1. Connection coefficients with the canonical sequence. On the grounds of (j2.ip , by per- 
forming straightforward computations we deduce that the coefficients c n ^ v (a) defined through 

n 

p n (x;a) = £ c n y(a)x v 

v=0 
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fulfill the recurrence relation 

c„+i,v(«) = (v + ce) c n . v (a) - (2v + 2a- l)c„ jV -i(a) , < v < n, , n^O, (2.6) 
under the convention of c n ,-i(«) = c nj „ + i(a) = 0, while (|2.2|) yields 

(2a + l)c„, v -i(a + 1) = -c„ + i,v(a) + a 2 c, i;V (a) , 1 ^ v ^ n+ 1 ,n ^ 0. (2.7) 
Lemma 2.2. 777e sequence p n (x;a) is then given by 

Proof. Let us set c, h v{tt) = (— 2) v (a + l/2) v c n>v (a) in order to deduce an explicit expression 
for c„ iV (a). The "triangular" relation (]2 .6(1 ensures another "triangular" relation for the new 
set of coefficients {c„ iV (a)}o^voi 

( c n+hv (a) =c„ jV _i(a)+ (v + a) 2 c„ jV (a) , ^ v ^n, ,0 0, 
[ c, !;0 (a) = a 2 " , c n . v (a) = , v ^n+l , n^O. 

The particular choice of n = 0, 1 or 2 in (j2.9j) furnishes the identities 

ci, (a) = a 2 , cij(a) = 1, c 2)0 (a) = a 4 , c 2! i(a) = 1 +2a(a+ 1), c 2 , 2 (a) = 1 , 
which show the validity of the identity 

for at least n = 0, 1 and 2. By finite induction, we will show that (j2.10p actually holds for any 
n e N. Indeed, according to (12. 9|) , it follows 



2(v + a) 2 ^ ^vV-l)^ +v r(2a + At )^, , ^ 2n+l 



+ ~ r(2« +At + v + l) (° + "> 



A 1 

v 



v! M % W r(2a + At + v + l) [a + N 
which corresponds to (j2.10p when n is replaced by n+ 1. □ 

The particular choice of a = gives a simpler expression for p n (x;0) than the one obtained 
in [16]. Moreover, in this case, the coefficients c„, v (0) coincide with the 0-modified Striling 
numbers of second kind, a particular case of the A-modified Stirling numbers of second kind 
treated in [7] - they correspond to the Jacobi-Stirling numbers expounded in [4] and explored 



6 ANA F. LOUREIRO, P. MARONI, AND S. YAKUBOVICH 

from a purely combinatorial point of view in [5]. Indeed, the set of numbers {c„ iV (a)} bear 
some resemblance with the Jacobi-Stirling numbers, triggering the problem of giving them some 
combinatorial significance. To avoid dispersion, we defer this study for a future work. 

2.2. Generating function. From this point forth we consider SR(a) > 0. Using relation 
(2.16.6.4) of [13] and the reciprocal formula (|1.5p we obtain the representation 

e"V = ^°°Tsinh(7rT)|r(a + /T)| 2 J S:, T (x)jT. (2.11) 

The absolute and uniform convergence of 

/ ' tW T 2 " +1 sinh(7TT)|r(« + /T)| 2 jT , m,«eN (2.12) 

JO ox m 

with respect to x ^xq > 0, easily verified by taking into account inequality 



d m K iT (x) 



dx n 



^e~ ST K m (xcosS), x>0, T>0, m G N (2.13) 



with 5 G (0, 7i/2), permits to interchange the order between the integral f)2 . 1 1 [> and the operator 
Hence, appealing to (jl.3p combined with (|2. . we derive 

*T = 7r3/2r(a + 1/2) I r 2n+l sinh(KT)\T(a + ir)\ 2 K it (x)dT (2.14) 

Therefore, recalling (jl.6p , the expression for polynomials p n (x;a), we obtain the integral rep- 
resentation 

p B (je;a) = (-1)" 7r3/2r(a + 1/2) y T 2 " +1 sinh(^)|r(a + /T)| 2 ^ T (x)jT , n G N„, (2.15) 

where the latter integral is absolutely convergent for all x > (see (12.121) ). 
By setting 

F a (u,x) = ^75— — e x x- a / cos(tw) t sinh(TTT) \F(a + ir)\ 2 K ir (x)dT (2.16) 

7r J ' z l (a + 1/2) jo 

and by taking into account that — ' ,n G No, is uniformly convergent by u in R via the 
absolutely convergent integral (12.151) , we obtain 

g J_a[U,X ) = Z ^)*X / T 2, 1+ l sinh(7rT)|r(a + /T)| 2^ TWjT = A! ^ ;a);?jGNo) 

«->o du zn n i i l T[a + 1/2) jo 

(2.17) 



7T 3 /2r(a + l/2) 

whereas 



M ->o du 2n+l K ' 



As a consequence, the series representation is allowed 

p„(x;a) 2n 
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We stress that 

^I^A = -e ( x- a ^- x x a F a (u,x), 
au L 

which, owing to (j2.4p , corresponds to a second order partial differential equation. Besides, the 
insertion of (j2. 19B into this latter relation yields (12. ip . Moreover, we have 

d 2n F a (u,x) 



,2n 



(-l) n e x x- a ^ n e- x x a F a (u,x) ,n£% 



2.3. Connection with the generalized Euler numbers. We integrate through (j2.15p multi- 
plied by e.~ 2x x a+e ~ l with respect to x over IR + , changing the order of integration in its right-hand 
side due to Fubini's theorem (use inequality (j2.13p ). Now, we invoke that 

: |r( £ + /T)| 2 



/ x e ~ X Q~ x K H {x)dx = 2~ e ^— 

Jo r 



(£+1/2) 



to obtain 

lim 

e->OJo 



/■CO 

lim / e~ 2x x a+e ~ l p n (x;a)dx 

e->OJo 



e^o 7rr(a + l/2)r(e + l/2) 
which yields 



poo I 1 \nni — U poo 

/ e- 2x x a - l p n (x;a)dx= i. / — / T 2 "\T(oc + iT)\ 2 dT , n G No- (2.20) 

jo v^r(a + 1/2) jo 

We recall the Generalized Euler polynomials E 2a (x) [fl Vol. III][SJ Q21 [15] defined through 

The particular choice of jc = a brings the identity 

(cosh( ? /2))- 2a = ££f(a)^. (2.22) 

On the other hand, [13 (1.104)] 

920;— 1 /• co 

(cosh( M /2))- 2a = — — — / cos(T M )|r(a + /T)| 2 ^T 
k 1 (2a) jo 

providing similarly to (I2.17P the equality 

yin / i \n j2a-\ poo 

api-w j )r ! * = i T2 " |r(a + ■ " e No - 

Hence, ([Z20]) together with (gjgp imply 

ja—l poo 

E^(oc) = =^(-1)" / e- 1 ^- 1 Pn (x;a)dx , R e N , (2.23) 
1(a) Jo 
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which, according to (|1.6j) . amounts to the same as 



E%(a) 



a-l nca „—X 



T a) 



sTe~ x i^dx , n G N . 

x 



Aware of the many formulas that can be found in the literature, we give here a formula that, 
as far as we are concerned, is new in the theory. 

Lemma 2.3. For any complex parameter a such that 91(a) > 0, the Generalized Euler poly- 
nomials of parameter 2a, E% a (x) are given by 



n- 1 



4=0 



n-k 



El a (x)=x» + Y j l £ (-2)- T (2a) T S(n-*,T) K 



(2.24) 



T=l 



where S(£,t) represent the Stirling numbers of second kind. 



Proof. The nth order derivative of the left-hand side of (j2.21h evaluated at the point t = 
furnishes an expression for E} x a (x). The Leibniz rule for the derivative of the product of two 
functions permits to write 



d" 



— — I e a 
dt n V e' + 1 I 



= 1 

f=0 v =° 



dt n-v \ e ' + 1 



2a 



neN, 







t=0 



According to the Faa di Bruno's formula for the kth order of derivative of the composition of 
functions (we refer to [3j Chapter III] for notation and a compendium of results for the Bell 
polynomials and the Faa di Bruno's formula), for any k ^ 1, we have 

2a k £ r t Q t Q t 

= E(a)M(- 1 ) M **,M(2'2'2'2'" 

=0 /* =1 r=0 
where B n ±{x\ ,xz, ■■■ ,x n -k+i), 1 ^k^n, represents the Bell polynomials evaluated at then-tuple 
(x\,X2, ■ ■ ■ ,x n ). On the basis of the properties of the Bell polynomials it follows (see [3] p. 135]) 



dt k 



e f + l 



(2.25) 



2' 2' 2' 2'' 



5^(1,1,1,1,1,...) 



S(k,n) 



with S(k,fi) representing the Stirling numbers of second kind [3] Chapter V]. The insertion of 
this information into (j2.25p , ensures 

2a k 

= £(a) M (-2)-"S(* J M) 

f=0 M=l 



A* V e' + 1 
whence the result. 

With x taken equal to a in (f2~2"4|) , it follows that £|" +1 (a) = 0, while 



□ 



4=1 



/< 

I 

T=l 



£ 2 2 « ( «) = + a"-* £ (-2)- T (2a) T S(*, t) , n e N , 
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are polynomials of degree n in a. We list a few examples of the (polynomial) sequence 

{E 2 "(a)} n ^ - 

Ef«(a) = -f £ 2a (a) = ^ + f 

plates _ 15a 3 _ 15a 2 _ a p2ce/„A _ 105a 4 i 105a 3 , 147a 2 , 17a 

c 6 l u l~ 8 8 2 8 v / 16 8 16 8 

2.4. Connection with the Bernoulli polynomials when a is a positive integer. Consider 
the (modified) falling factorial of a complex number x, 



' m— 1 

n 

* (7=0 



1(x-a 2 )] , neN, and [x] := 1. 



This is actually a polynomial of degree n whose coefficients in the canonical basis are the 
so-called 0-modified Stirling numbers of first kind, denoted by so(n,v), and treated in [7J. 
Regarding their importance for the subsequent developments, we recall 

[x} n = J £s (n,v)x v , «6N . (2.26) 

v=0 

This gives grounds for writing |T(a + /t)| 2 , whenever a = m G N, as follows: 
As a consequence, the generating function F m (u,x) given by (I2.16P becomes 

m poo 

F m (u,x) = C m % x- m e x Y (-l) m+a s (m, a) / t 2a cos(zu)K iT (x) d% 
(7~o Jo 

which amounts to the same as 



m _2 ~\2a , . 

F m (t,x)=C m (-1)'" x- m £ 1L ? (m,o) ^(e- 2xsnih ( " /2) ) (2.27) 



because of the equality [TBI (112)] 

J t 2 " cos(m)K H (x) dx = y > 2 — & -*™Mu) j „ G Nq _ 

Concomitantly, the polynomials p n (x;m) for wi £ N are related to p„(x;0) through 

7T 2 m 

x> ;i (x;m) =C m (-l) m y £ Sb(m, o)p n+a (x;Q) , n > 0. (2.28) 



(7=0 
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2 1 

with C m = — • — ; — ;— ; — ■ Recalling [161 (3.14)1, it follows a connection with the Bernoulli 

n 2'»(l/2) m ;J 

polynomials evaluated at 1 



2 

2 2(«+l) e * j^^j /1-/T N 

'0^«,o;— 

(T=0 



71 m — 2 2( " +1 V f°° Kt(x) (\ — iz\ 

x m p n {x;m)=C m (-1)'"— V T (m,a) / T— B 2n+ i — ) dx , O 0. 

2 ^ (2n + 1)7T/ Jo x \ 2 J 

(2.29) 

Such equality has the merit of deriving analog properties for p n {x;m), when m G N, to those 
obtained in [15]. Namely, the connection with the Bernoulli or Euler numbers as well as the 
Bernoulli polynomials. For further considerations, we refer to [16, (3.3), (3.13)-(3.14) and 
(3.17)]. As a consequence of the aforementioned, we have the following result, which, as far 
as we are concerned, is new in the theory. 

Lemma 2.4. The Generalized Euler numbers of parameter 2m and the Bernoulli numbers are 
related by 

£ " (m) = (£.-1)1 £^-?o(",")«^ , n e N„. 

Proof. The relation (3.3) pointed out in [15] furnishes an identity between the Bernoulli num- 
bers and the polynomials p n (x;0) by means of an integral 

dx 



n r°° dx 
B 2n = x _ 22n j Q e- 2x Pn(x;0)— , n € Ni 



The combination of this latter with (j2.23h together with (j2.28H brings to light the desired 
equality. □ 

The forthcoming developments leave aside aspects of the polynomial sequence {p n {x;a)} n ^>o 
itself in order to embrace outer consequences. The study will be essentially focused on the dual 
sequence corresponding to the monic polynomial sequence, hereafter {P n (x;a)} n ^o, obtained 
from {p,i{x;a)} n ^o by dividing each of its elements by the respective leading coefficient: 

P n (x; a) = . , ~r~ — — — p n (x; a) = — 1 - j- e v x~ a s^ n e~ x x a , n G N . (2.30) 
y } (-2)" (a + 1/2),/ v ; 2"(a + l/2)„ v ' 



3. The monic polynomial sequence and its corresponding dual sequence 

As, from now on, we will deal exclusively with the monic polynomial sequences {P n (x;a)}„^o, 
we recreate the needed properties based on those already obtained for the original sequence. 
Indeed, in accordance with (12. lj) , {P„(x; a)}„^o fulfills 

x 2 P^(x\a) -x(lx- \ -2a\ P' n (x;a)- ((2a + l)x- a 2 ) P„(x;a) = -(2n + 2a + l)P n+l (x;a) 

(3.1) 
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for all n £ No, whereas the equations (I2.2j) - (j2.5j) yield 

a 2 

x P n (x;a + l)=P n+l (x;a) + —— — — P n (x;a), n e No. (3.2) 

(2n + 2a + lj 

In the next §3. II we will ferret out properties of the corresponding dual sequence, {u n (a)} n ^>o, 
which will trigger some interesting results, namely, the positive-definite character of the canon- 
ical element uq(oc) as it will be shown in 



3.1. The dual sequence. The differential and difference equations satisfied by the MPS 
{P n (x;cc)} n ^o enable the differential and difference equations fulfilled by the corresponding 
dual sequence {u n (a)} n ^o- 

Lemma 3.1. The dual sequence {u n (a)} n ^o associated to {P n (x;a)} n ^o is such that 

((x 2 w (a))' + x(2x-(l+2a))w (a)) - ((1 +2a)x- a 2 )u (a) = (3.3) 

{(x 2 u n+ \(a))' + x{2x- (1 +2a))w„ + i(a)^ - ((1 + 2a)x- a 2 )w, 1+ i(a) = (2n + 2a + \)u n (a) 

(3.4) 

for n G No- Moreover the moments of the canonical form uq are given by 

Proof. The action of uq over the relation (|3.1j) furnishes 

{u (a),-x 2 P"(x)+x(2x-l-2aJ P' n (x) + ((2a + l)x- a 2 ) P n (x))=0, O 0, 
which, by transposition, in accordance with (jl.8|) , becomes 

(-(x 2 u (a))" - (x{2x- I -2a) u (a))' + ^(2a + l)x- a 2 ^w (a) , P„(x)^=0, n^O, 
whence (j3.3p . Likewise, the action of u^+i over (13. ip provides 

u k+ i(a) , -x 2 P"{x)+x(2x-\-2a s j P,'(x) + {{2a + l)x- a 2 } P„(x)\ = (2n+2a+l) S„ jk , 0, 
and, again, by transposition, this latter leads to 

(-(x 2 u k+ \(a))" - (x(2x- 1 -2a) Uk+i{a))' + ((2a + l)x- a 2 ) u k+ \ (a) , P„{x)) = (2n + 2a + 1) 8 n ^-, n 
therefore, according to fll.7f> . we derive fl3.4p .On the other hand, the action of both sides of 



(I3.3I) over each element of the sequence {x"}„^o leads to the following difference equation 
having the moments of uo(a) as solution 

(2n + 2a + \)(u Q (a)) n+l = -(n + a) 2 (u (a)) n , n e No. 

providing Q.5p . □ 
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In addition, according to (II. 7j) , we have as well 

xu (a) = Y,( x u (a),P n (x;a + l))u n (a + l). 

Considering (II. 8p and (j3.2[) . it follows 



2 2 

(x u (a),P„(x;a + l)) = (« («) , P„+i (*; a) + — - 2a - - P„ (x; a)) = 2n + 2a + l 5nSh 
that holds for any n S No. Therefore the two forms uo(a) and uo(a + 1) are related by 

cc 

x u Q {a) = 2a + l «o(« + l) • (3.6) 

The question on whether the MPS {P„(x; a)}„^o can be orthogonal arises in a natural way and, 
concomitantly, on whether the form uo(a) is or is not regular. 

As a matter of fact, we recall that a form v is said to be regular if we can associate with it a 
PS {<2„}n>o such that (v, Q n Q m ) = k n 8 n>m with k n ^0 for all n,m eN 0HQ|. The PS {gnWo 
is then said to be orthogonal with respect to v and we can assume the system (of orthogonal 
polynomials) to be monic. Therefore, there exists a dual sequence {v„}„^o and the original 
form is proportional to vrj. Furthermore, it holds 

v n +i = ((vo,GLi(-))) _1 G»+i(*)vo (3.7) 
Moreover, when v G is regular, let <I> be a polynomial such that <£v = 0, then <I> = [91 [10] . 

In this case we call this unique MPS {Q n {x)}n^o as monic orthogonal polynomial sequence - 
hereafter MOPS - and it can be characterized by the popular second order recurrence relation 

{Qq(x;o) = \ ; Q 1 (x;cc)=x-po(a) 

\Q n+2 (x;a) = (x-p n+1 (a))Q n+ i(x;a)-Y n+1 {a)Q n (x;a) , n£N Q . 

Here, we consider the dependence on the complex parameter a to avoid repetition of similar 
formulas. We will systematically refer to the pair (/3„(a),y„ + i(a))„^o as the recurrence coeffi- 
cients of {Qn{x;cc)} n ^o, necessarily fulfilling y n+ \{a) / 0, n ^ 0. See [51 [TU] for notation and 
a compendium of results about algebraic properties of orthogonal polynomial sequences along 
with regular forms. 

Somehow expected, we have the following result. 

Lemma 3.2. The MPS {P n (x;a)} n ^o cannot be regularly orthogonal. 

Proof. Under the assumption that {P,,(x; ce)}„^o defined through (j2.30p is orthogonal, we may 
insert relation p. 71) , with v n replaced by u n (a) and Q n by P n (-;a), in (j3.4p which provides 



(x 2 P n+ \(x;a)uo(cc)Y +x(2x- (1 +2a))P n+ i(x;a)u (a) 

-((1 + 2a)x- a 2 )P n+ i(x;a)u (a) = X n (a)P n (x;a)uo(a), n e No, 

(3-9) 
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where X„(a) = (2n + 2a + 1) < uo(a),P 2 +l (x;a) > (< uo(a),P 2 (x; a) >) , for all n e No. 
Taking into account (13. 3p , then (13. 9 p implies the differential equation 

2P' 1+i (x; a)(x 2 u (a))' + (x 2 P" +l (x;a) +x(2x-(l + 2a)P' l+l (x; cm u (a) = X n (a)P n (x; a)u (a), n 
When the equation obtained by the particular choice of n = 0, i.e. 

2(x 2 uo(a))' + x(2x - (1 + 2a)w (a) = Ao(a)wo(a), 
is inserted into the original equation, we obtain 

(x 2 P" +l (x;a)^u (a) = {x n (a)P n {x\a) - h ) (a)P l n+l {x\a) s ju Ql (a). 
The regularity of uo(a) would now imply the condition 

x 2 P^ +l (x;a) = X n (a)P n (x;a)-h (a)P / „ +1 (x;a) , «>0, 

contradicting degP n (x; a) = n, and therefore crumbling the possibility of {P n (x;a)} to be or- 
thogonal. □ 

Despite the non-(regular)orthogonality of {P„(x; cc)}„^o with respect to the form uq(cc), we 
cannot exclude the existence of an orthogonal polynomial sequence, say {Q n (x; cc)} n ^o, with 
respect to uo(a), which amounts to the same as ensuring the regularity of uo(a). This question 
is handled in the next section. 

3.2. About the regularity of uq. We begin by rewriting (13. 3p as follows 

(0M O (a)) / + V^Ko(a)) +£wo(a) = (3.10) 
with 

(j>(x) := (j)(x;a) =x 2 , y(x) := y(x;a) =x(2x-2a - 1), %(x) := x{ x '^ a ) = -(2a + l)x + a 2 . 

(3.H) 

Actually, while seeking an integral representation for uo(a), we realize that it is indeed regular. 

Lemma 3.3. For any positive real value of the parameter a, the form uo(a) is positive definite 
(therefore regular) admitting the integral representation 

(u ,f(x)) = 7 'J / f(x) e- x x a - l K (x) dx , V/ € ^. (3.12) 

yjn 1 {a) z Jo 

Proof. We seek a function U(x) :=U(x;a) such that (j3. 12H holds in a certain domain C. Since 
(uq, 1) = 1 0, we must have 



J U{x)dx= 1^0 (3.13) 
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By virtue of ()3.10p , we have, for any / G 3? 

= (({<$>{x)u )' + y(x)u^ + x{x)u Q J{x)) = (u ,$(x)f"(x) + ¥ (x)f'(x) + x(x)f(x)) 

{((f>(x)U(x))" + (Y(x)U(x))' + %(x)U(x)) f{x)dx 
- (4>{x)U(x)f(x) - ((j>(x)U(x))'f(x) - Y(x)U(x)f(xj) \ c 
therefore, U(x) is a function simultaneously fulfilling 

J c {($(x)U(x))" + ( ¥ (x)U(x)y + X (x)U(x))f(x)dx = , Vfe& (3.14) 

((t>(x)U(x)f(x)-($(x)U(x)yf(x)- V (x)U(x)f(x))\ c = , Vfe&. (3.15) 
The first equation implies 

(4>(x)U(x))" + Mx)U(x))'+x(x)U(x) = Xg(x) 

where X is a complex number and g{x) ^ is a function representing the null form, that is, a 
function such that 

Jg(x)f(x)dx = , Vfe&. 
We begin by choosing A = and we search a regular solution of the differential equation 

(0 (x)U(x))" + (y(x)U(x))' + x(x)U(x) = 
Upon the change of variable U(x) = e~ x x a ~ l y(x) we have 



x 2 y" (x)+xy'(x)-x 2 y(x)=0 d (y (x) \ 







whose general solution is: y(x) = c\Iq(x) + C2Kq{x) , x ^ 0, for some arbitrary constants c\,C2 
and y(x) = when x < FjQ. As a consequence U(x) = e _A jc a_1 jci/o(*) +C2^oW| , x ^ 0. 
Insofar as U(x) must be a rapidly decreasing sequence (that is, such that lim f(x)U(x) = 

for any polynomial /), we set ci = and C2 / and we write C7 (x) = C2&~ x x a Kq{x) , x £C = 
]0, +°°[ with C2 set in order to realize (I3.13p , which amounts to the same as C2 = 2 J^Y^Jy ■ ' n 

this case U(x) also fulfills the condition (j3.13p as long as a > because e~ x x a ~ 1 Ko(x)dx 
is a stricktly positive convergent integral. Besides, since K v (x) has the asymptotic behaviour 
with respect to x [TJ Vol. 1 1] [19] 

/ JT \ 1/2 

X v (jc)=(-) e-*[l +0(1/*)], x^+°o 



>2x/ 

^ v (jc) = <9(x~ 9?(v) ) ,Xb(jc) = 0(-logJc) , x^0, 

then we immediately deduce that (j3.15p (13 . 15[) is fulfilled by any element of the PS {x"} n ^,o 
(that spans and hence it is fulfilled by any element / of & . Moreover, for every polyno- 
mial g(x) that is not identically zero and is non-negative for all real x we have (uo,g(x)) = 
2 a F(a + 1/2) r+°° 

— — - — — — / g(x) e x x a Ko{x)dx>Q and therefore uq is positive-definite form. This 

V k r(a) 2 Jo 
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implies that uq has real moments and a corresponding MOPS exists (i.e., uq is a regular 
form). □ 

The regularity of uq((x) raises the problem of characterizing a corresponding orthogonal poly- 
nomial sequence, say {Q n (x; oc)} n ^o, whenever a is a positive real number. From this point 
forth we consider a € IR + . 

Entailed in the question of characterizing the arisen MOPS {Q n (x;a)} n -^o (with respect to 
wrj(a)), comes out the problem of determining the associated recurrence coefficients or even 
a differential equation fulfilled by {Q n (x; ot)}„^o or an analogue of Rodrigues' formula or a 
generating function. However, this has revealed to be a tricky problem to solve, reminding 
other open problems such as the one posed by Prudnikov [14]. Indeed, the difficult part of 
this problem is connected to the fact that the regular form uo(a) is solution of a differential 
equation of order higher than one. As far as we can tell the gap in the theory concerned to 
this problem, avoids us to attain further results. Notwithstanding this, based on the work [11J, 
it is possible to reach a finite-type relation (of order two actually) between the two MOPSs 
{Q n {x; a)} n j»o and {Q n (x;a + 1)}„>o- Again, this adds very few to the answer of the problem. 

The recurrence coefficients (P„(a),y n+ i(a))„^o associated to the second order recurrence re- 
lation (|3.8p fulfilled by {Q n (-',(x)}n^o mav be successively computed either by means of the 
weight function or using the Hankel determinant of the moments of uo(a) [2]. Making use of 
the first process, for each positive integer n we have 



(u (a),xQl(x;a)) 



f + °° 

/ Q n (x;a) e- x x a K (x) dx 
Jo 



= v y 7::^r: = (3.16) 



(uo(a),Ql(x;a)) 



Q 2 n (x;a) e - V _1 Xb(jc) dx 

r+o<> 



r+<*> 

_ («,(«),&,(*<»)) _ I <=-*°-'*i>M <•* 

Jo 

According to this latter, we list the first elements 

RJrv\ - _«L ft, <n\ - »(2«(«+4)+7)(o:(2o!+5)+4)+4 

POl" J — 20 ;+l ' Pl\ u ) — (2a+l)(2a+5)(2a(a+2)+l) 

/ x _ a 2 (2o;(«+2)+l) / x _ 4(«+l) 2 (2a+l)(a(a(2o;(ct(2o:(a+12)+l 13)+262)+613)+325)+51) 

(2a+l) 2 (2a+3) ' ' 2( ^ U ' 1 (2a+3)(2a+5) 2 (2a+7)(2a(a+2)+l) 2 ' 

The determination of recurrence coefficients of higher order has indeed revealed to be a ticklish 
problem. 
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